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Abstract 

The relativistic kinetic equations (RKE) for lepton plasma in the presence of a strong 
external magnetic field are derived in Vlasov approximation. The new RKE for the elec- 
tron spin distribution function includes the weak interaction with neutrinos originated by 
the axial vector current (~ ca) and provided by the parity nonconservation. In a polarized 
electron gas Bloch equation describing the evolution of the magnetization density pertur- 
bation is derived from the electron spin RKE being modified in the presence of neutrino 
fluxes. Such modified hydrodynamical equation allows to obtain the new dispersion equa- 
tion in a magnetized plasma from which the neutrino driven instability of spin waves can 
be found. It is shown that this instability is more efficient e.g. in a magnetized supernova 
than the analogous one for Langmuir waves enhanced in an isotropic plasma. 
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1 Introduction 



When he conjectured the existence of the neutrino, PauH imposed very stringent bounds 
on its electrical neutrality. 

Nevertheless, the direct (=weak) interaction of neutrino with electrons and positrons 
shows a principle possibility of the interaction of neutrino with the electromagnetic field. 
In vacuum, however, this interaction is negligible. The situation is extremely changed in 
media with free carriers of the electric charge such as a dense plasma of metals, stars, the 
lepton plasma of the early Universe, etc. 

It was shown in |^, ^ that the electrodynamics of neutrinos is changed in such 
media comparing with the electrodynamics of neutrinos in vacuum in such a way that the 
contribution of the neutrino electromagnetic vertex in medium cannot be regarded as a 
small correction like the case of radiative corrections in vacuum. The main distinction 
is the appearance of the induced electric charge of neutrino which is proportional to the 
density of free carriers of electricity Namely this leads to the new qualitative effect: 
the appearance of the long-distance forces inevitably leads to the collective interactions 
of neutrino with an ensemble of charged particles via electromagnetic field. 

The approximation of the neutrino propagation in an isotropic plasma considered in 
dl, 1^, § and then in Q has some natural bounds for applications in astrophysics. 
Therefore one needs to have the self-consistent system of kinetic equations of the most 
general kind applicable for other astrophysical plasmas including magnetized stars. 

Really this turns out to be a new problem unsolved before in the neutrino kinetics 
while in the one-particle Schrodinger equation approach one can show the importance of 
the axial vector potential for the neutrino propagation in a magnetized plasma (Va ~< 
■0e7i75'0e >~ Bi) that provides the interaction of neutrino with the magnetic field B 
without any neutrino magnetic moment Q and could explain, e.g. the neutron star kick 

i- 

In the present work we take into account the spin interactions of electrons that are 
important for neutrino propagation in a dense degenerate electron gas polarized by the 
external magnetic field. 

In subsection 2.1. we discuss the equilibrium state of an electron gas polarized by the 
external magnetic field Bq ^ 0. 

Then in subsection 2.2 we present the full set of the coupled Relativistic Kinetic Equa- 
tions (RKE) for electrons and neutrinos in the lepton plasma with an external magnetic 
field including the self-consistent electromagnetic field contribution. These RKE's are 
derived using the Bogolyubov method analogously to the way described in details for an 
isotropic plasma in |^]. 

In subsection 2.3 we complete derivation of RKE's recasting them for the gauge in- 
variant distribution functions and prove the lepton current conservation, 
djj!^\x,t)/dxf^ = 0, in a magnetized plasma both for neutrinos and electrons. 

We conclude this section deriving in subsection 2.4 the dispersion equation for the 
lepton density perturbations in an isotropic plasma that coincides with results [0, 10|. 



In section 3 we generalize the Bloch equation for the evolution of the magnetization 
density perturbation in the presence of neutrino fluxes. In section 4, neglecting spatial 
dispersion we obtain the increment of spin waves excited by the neutrino beam and com- 



^ There appears also induced magnetic moment of neutrino cross-sections of the neutrino scattering off 
nuclei are modified in plasmas especially in low energy region IpTand so on. 
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pare it with the case of the neutrino driven instabihty of plasma waves in an isotropic 
plasma. 

In section 5 we discuss results and give some estimates of relevant parameters in a 
polarized medium. 

2 Relativistic Kinetic Equations for lepton plasma 
in Vlasov approximation 

The kinetic equations in the Standard Model (SM) are derived from the quantum Liouville 
equation for the non-equilibrium statistical operator p{t) using the Bogolyubov method 
with the total interaction Hamiltonian given by the Feynman diagrams for the neutrino 
scattering off electrons and the usual ee-scattering in QED Q. For a magnetized plasma 
we do not consider i^i/-scattering and neglect also weak ee-scattering comparing with the 
electromagnetic interaction of charged particles through the photon exchange. 

The one-particle density matrix /X''(p,x,t) = Ek e'''''rr(p(t)6pi",^/2,r^p+k/2,r') 
termined as the statistical average of the product of creation (6pi"^) and annihilation 
{b^pi^r') operators for the lepton a = e,^. In the case of electrons it takes the form 

/S(pe,X,t) =/(^)(p„X,t)^+5f)(p.,X,t)M^ , 

where A = ±1 is the eigenvalue of the spin projection on an external magnetic field. 
E.g. for the magnetic field B = (0,0, i?o) only the spin component T,z commutes with 
the Hamiltonian for electrons, [S^, ?i(^)]=0, hence ((T^)a'a = ^^\'\- In the absence of 
neutrino magnetic moment the one-particle density matrix for massless Dirac neutrinos 
corresponds to the active (left-handed, r = — 1) neutrino distribution function only, 

/ir^q,x,t) = il::^/H(q,x,t), 

where in an uniform medium the helicity r = ±1 conserves, {(Jiqi)r'r = qTS^'r- 

The number density distribution functions are related with the lepton densities n^(x, t) = 
/(dV(27r)^)/(^Hp,x,t) for electrons and n^(x,t) = /(d^g/(27r)3)/(^)(q, x, i) for neutri- 
nos. 

2.1 Equilibrium state in a polarized electron gas 

In the linear approximation for a slightly inhomogeneous medium the distribution func- 
tions have the form 

/M(q,x,t)=/('^)(q)+5/H(q,x,t) , 
/(^)(p,x,t)=/('=)(6p)+<5/(-)(p,x,t) , 

S(^)(p,x,t) = s(^)(ep)+5S('=)(p,x,t) , (1) 

where the neutrino background is given by /o'^^(q) not being in equilibrium with plasma 
environment (e.g. for the powerful neutrino flux outside of the SN neutrinosphere) , and we 
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consider the uniform polarized equilibrium electron gas for which the one-particle Wigner 
density matrix takes the form 

&\e,) = '-ffi^\e,)+^^sf\e,). (2) 

Here fW{ep)=ge [exp(ep — C)/T) + 1] ^ is the equilibrium Fermi function; = 2 is the 
Lande factor; 



(3) 



is the equilibrium spin distribution; 7 = Ep/rrie is the electron gamma factor, and we use 
hereafter units h = c = 1. 

Note that in @ we assumed the quasi-classical electron spectrum in a realistic external 
magnetic field Bq obeying the inequality ei?o ^ T^, when the Landau spectrum 

e{X,p„n) = {ml+pl+ \ e \ 5o(2n + 1 - A))!/^ (4) 

reduces to 

tw A I e I i?o 

= ep — , (5) 

where A = ±1 and we changed | e | Bo[2n + 1) = for large Landau numbers n ^ 1. 
Thus, the full spectrum (|5|) contains the small paramagnetic (spin) correction to the 
continuous spectrum Sp = \Jp^ + ra\ . 

One can easily check that in this quasi-classical limit the exact expression for magne- 
tization of the electron gas, 

Mf^ = \fiB\< Aljl5i^ >o=| /iB I E dp,Tr[a,P'\e{X,p,,n))] = 



n=0 



=1 /^B I £ * ^!'f\2^' I dp^Y.^fx\\^i^^Pz,n)) , (6) 

takes the quasi-classical form 

Mf =1 MB I / ^,sf{'v) = -2mI^o, / D{ep)^^^dep , (7) 

where the spin distribution S^p {ep) is given by eq. d) and D{ep) = pmeli^-nf. In 
particular, in the non-relativistic (NR) limit this background magnetization corresponds 
to the spin paramagnetism of the free electron gas in metal 1 11 1 . 

Note that for the degenerate electron gas from the quasi-classical Eq. (^) one obtains 
the same value Uf^ =1 /is II e I Bqpf^/2t:'^ as it follows from the exact quantum eq. (g) 
for an arbitrary strong magnetic field for which electrons populating the main Landau 
level (n=0) contribute only. 

Note also that this magnetization determines the axial vector potential of a probing 
neutrino in a magnetized plasma too |§, 



Va = -GfV2ca'^^ 



(0) 



Q I f^B 
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that changes the spectrum of the ultrarelativistic neutrino, £q = q + V + Va, comparing 
with the standard one in an isotropic medium (e.g. in the Sun), Eg = q + V , and modifies 
the neutrino oscihations in a magnetized SN described by the one-particle Schrodinger 
equation for two neutrino species [S, 0], 



V + Va- C25 
S2S 



S2S 




(8) 



where C2 = cos 20, S2 = sin 20, 6 = Arn^/AE are the neutrino mixing parameters; Vb is 
the sterile neutrino wave function. 

In our kinetic approach for the SM lepton plasma we do not consider oscillations 
(s2 = 0) assuming the massless spectrum £q = q for active neutrino of the given kind Va-, 
a = e, fj,, T. 



2.2 The master kinetic equations for perturbations in plasma 

Within the linear approximation (|^) the RKE for the Lorentz-invariant Wigner number 
density distribution functions Sf^"'\p,x,t) = Tr 5/^ \p,x, t) 



take the covariant forms 



dxa J (27r)3 dxj dqj 



GpCAmeq^ 



(2vr)= 

d^Pe d5a\,'\pe,^,t)df^''\q) 



+ 



V2 J i2^?ep. 

for neutrinos and 

a(^/(^)(p,x,t) 



dxj 



dqj 



, 



(9) 



dx. 



+ e6Fjf,{^,t)p^ 



' iu P 



dpj dpj 

+GfcvV2 f {p^q,)d6f^'^\q,^,t)dfS,'\ep) ^ 

J (27r)3 £g dxj dpj 

J [zir)'^ Eg opj oxj 



+ 



(10) 



for electrons where cy = 2^ it 0.5, ca = =f0.5 are the vector and axial weak couplings with 
upper (lower) sign for electron (muon or tau) neutrinos correspondingly, ^ = sin^ 9w is 
the Weinberg parameter; the tensor Fj'^ = Cj^iBQi corresponds to the external magnetic 
field term. 

Here we have introduced in eq. (^) the 4-vector a}f^(p,x, i) = a^p)S'*(p, x, t) that is 
the statistical generalization Q of the Pauli-Lubahski 4-vector [jl4 



a^ip) = (ar(p)O) 
and has the components 



pC.^-, p(pC) 

me' me{Ep + me) 



ps^'Hp,^,t) (,) p(psW(p,x,t)) 

,a (^p,x,rj-i 1 — 

me{Ee + me) 



rup 



(11) 
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Note that in the presence of the external magnetic field Bq the neutrino parts of these 
Vlasov equations differ from the analogous ones used in due to the inclusion of axial 
vector interactions (~ ca) in the total i^e-scattering amplitude that allows us to account 
for collective interactions of and Vr -neutrinos too for which the vector coupling cy is 
smah, cv = 2i- 0.5 ^ 0. 

To check general symmetries we temporarily refuse the linearization Eq. ([l|) and write 
down the RKE for the total electron spin distribution function S'-'^'*(p, x, t) accounting for 
the weak interaction of electrons with neutrinos: 



a5f)(p,x,t) asf(p,x,t) , p^^dgHp^^ 

+ V 1- el'j^[-x., t) h 



dt 



5x 



dpj 



i?,(x, t)(vS(^) (p, X, t)) - ^;,(E( x,t)S(^)(p,x,t)) ^ [S('^)(p,x,t) x B(x,t)] 



rup 

H ^iki 



1 + 7 

A (^X(p))/H(q,x,t)5;^)(p,x,t) + 
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+ 



+ 



+ 



CATUe 



(2vr)3 sq 
d\ (P)) (p, X, t) (q, X, t) 



{2-Kf 



dpj 



dxj 



+ 



GfV2cv f d\ (p^(7^)55f)(p,x,t)5/H(q,x,t 



(27r) 



, GfV2cv , 



Eq dpj dXj 

d?q PrnQk ( (9/(''Hq,x,t)^ 



+ 



(2vr)3 sq 



dx. 



d^q g^pP<,(p) (q, X, t) 

(27r)3 Eq dXm 



eijnSf\p,:>ii,t) = . 



:i2) 



This is the relativistic generalization in SM of the kinetic equation for spin waves in 
nonferromagnetic metals |13|. But it should be kept in mind that we consider here the 
Fermi gas of free electrons in contrast with the quasi-particle approach for metals and 
also neglected exchange interactions (neglecting exchange Feynman diagrams for the ee- 
scattering means here that the long-range forces are dominant). 

The system of RKE's is completed by the Maxwell equations for the electromagnetic 
field F^y{x, t) = 5F^^{x, t) + F^^J accounting for the spin wave contribution, 



iaB(x,t) 
c dt 



[V xE(x,t)] , 



and 



[Vx[VxE(x,f)]]+ 2 



1 92E(x, t) in d 



+IJ-B 



d^p [V X S(")(p,x,t)] 
(27r)3 ^ 



■7^ e 



d^p 



v/(^)(p,x,i) + 



c2 at V J (27r) 
d'p ,^^^[yxS(^)(p,x,t)] 



(2vr) 



7 + 1 



(13) 



In eqs. (p|)-(p^ e = — | e | is the electric charge of the electron; Gi? is the Fermi constant; 
fj-B = eh/2meC is the Bohr magneton; v = p/e^ is the electron velocity; the latin indices 
run i,j = 1,2,3 and the greek ones are n = 0,1,2,3 for scalars written in the Feynman 
metrics, A^B^^ = AqBo - AiBi. 
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In the non-relativistic (NR) limit | v |^ 1 (7 — > 1) the last term drops out and eq. 
(13) coincides with the Maxwell equation in a magnetized medium: 



[V X (B - 47rM)] = 47rj + 



dB 
'dt 



(14) 



where M(x,t) =| \ /f^^pS(p,x, t)/(27r)'^ is the magnetization density of NR plasma. 

Multiplying eq. (|12|) by the energy Sp and the spin distribution function ^^(p, x, t) 
we obtain the covariant RKE for the Lorentz-invariant product a^(p, x, t)a^(p, x, t) = 
— (S('^)(p,x,t))^ where the 4-vector a|f^(p,x) is given by (11). 

Note that for the uniform electron beam {sl^\p,t) = Si{t){2-K)'^neo5^'^\p — po)) 
omitting the neutrino term and integrating the spin RKE (|l^) over d^p we obtain the 
one-particle spin evolution equation 



dt 



70 



70 + 1 



S X 



E X vq 



that turns out to be exactly the Bargman-Mishel-Telegdi equation for the electron spin 



motion with the normal magnetic moment jjb [14]. Here vq = p/^po ^he velocity of the 
electron beam, 70 = epg/nie is its 7-factor. 

In general, it is possible to generalize eq. (|^) adding the electromagnetic scattering 
of electrons through anomalous magnetic moment (Schwinger correction fj,' = {a/2'K)(jLB) 
that could lead to additional terms similar to electromagnetic terms in the neutrino spin 
RKE |1|]. 

Below we consider some important particular cases of the master equations checking 
their consistency with the known results 0, IC, 12 1. 



2.3 Lepton current conservation dj^°''{y.^t)/dx^ = 

Let us rewrite eq. (^ as the classical equation 

a5/M(q,x,i) , .a5/H(q,x,i) , .df^\ci) 



+ x- 



+ q- 



, 



dt 5x 3q 

where x = n = q/g is the velocity of massless neutrino, the derivative q is given by 



(15) 



q= + V 



GF\/2cy(5n(^)(x,t) + ^CAMo(x,t) 
v2 



-Gp 



V2cvV[n5i^'\^,t)] + ^V[n(^A(^)(x,t)] 



(16) 



Here 5jt^(x,t) = (5n(^)(x, f); Jj^*^) (x, t)) = j{d^p/{2TTf){p^/ep)5f^^\p,^,t) is the four- 
vector of the electron current density perturbation ; 

(5^}f'*(x, i) = me /(ii^p/(27r)^)(5a^'^^ (p, x, i)/ep is the axial four-vector of the spin density 
perturbation^; the four- vector spin distribution 5a|f^(p,x, t) is given by (0);V = di = 
d/dx^ = —9*. 



^In NR plasma the spin density is given by the 3-vector component only, (5A^(x, t) — (0, 6S{k, t)) 
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Analogously the electron RKE (|lO|) can be rewritten as the quasi-classical one, 

at ax \ ap / 

where the particle number distribution function is obtained via the summing over spin 
quantum numbers, 5/'''^^(p, x, t) = Tr{5f^^\'p, x, t); x = v = p/ffp is the electron velocity; 



^ = [ e(5E(x,f) + [v X B(x,t)]) + Gi.\/2cy[— V<5ji^)(x,t) 



+GfV2cA-^V5i(-)(x,t)a^(p)(afc),y (18) 

is the force matrix which accounts for the Lorentz force with electromagnetic fields 5E, 
B = Bq + (5B as well as the weak interaction terms ~ Gi;' 0; 

= ^-ff%e,) + (19) 

is the equilibrium density matrix as given in (^) with the ort b^'^^ = Bo/i?o separated 
from the isotropic distribution function S^^^\ep) = — (| /is | BQ/'y)df^^^^ /dSp given by 
i). 

Substituting ([l^), ( p!9[ ) into (|l^ one can easily check that the electron RKE (10) takes 
the form 

— dt — + " — K — + '^~p^^ + '^^-^ i; — — + 



dpj 

where we introduced the unit polarization four-vector a'^(p) = af{p)bf^\ a^a^^ = —1. 

On first glance both the neutrino ([l5|) and electron (17) RKE's do not obey the lepton 
number conservation law, djj^\:x.,t)/dx^ 7^ 0, (due to the last terms in (|l^ ) and the last 



three terms in (|20| ) correspondingly). 

This non-conservation is true for the gauge non-invariant Wigner distributions f^""^ (p, x, t) 
/('^)(p,x, t) = / e*Pyd'^y/'^''^(x — y/2, x + y/2, t) we used above, for which the distribution 
function in the coordinate representation (within integrand), 



/»(xi,X2,t) = Tr (p(t)v&(»)+(x2)^'(")(xi)) , 



^ Such tems can be also obtained from the weak interaction Hamiltonian for a probe electron moving in a 
neutrino medium 



^ /^"a(p)7''wv(p) (^) , ^ /t; "a(p)7''75Wa'(p) r -M/ ,n 
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is not invariant with respect to gauge transformation of the corresponding interaction 
fields. 

For example, in QED plasma the kinetic equation for the gauge non-invariant distri- 
bution of charged particles /('^^(xi, X2, t) derived from the quantum Liouville equation by 
the same Bogolyubov method does not obey the electric current conservation law since 
the force term depends on the electromagnetic potentials Afj_{x,t) which do not enter as 
combinations expressed via field strengths, E, B [^]. 

The recasting of such RKE for the gauge-invariant Wigner distribution 

/('=)(p,x,t) = /»(p + eA(x,t),x,t)= I dhe'Pyf^^\x-y/2,^ + y/2,t) , (21) 

where the gauge invariant distribution in the coordinate representation /(^)(xi,X2,t) is 
connected with the gauge non-invariant one, 

/('')(xi,X2,t) = Tr (p(t)l'(^)+(x2)^(^)(xi)), by the important phase factor [0]: 



/(^)(xi,X2,t) = exp 



ze(x2-xi)^ d,^A(^X2 + C(xi - X2),t 



/»(xi,X2,t) , (22) 



allows to obtain the usual Lorentz force term in the standard Boltzman equation for 
charged particles []l6|: 



ot ax op 

for which, of course, the electric current jl^\-x.,t) = J (i^p(p^/ep)/*^'^)(p, x, t)/(27r)^ is 
conserved, djj^^ /dx^ = 0. 

Such conservation is stipulated by the presence of the phase factor in the distribution 
(22) which in turn is invariant under the standard gauge transformation (with an arbitrary 
gauge function x(x, t) obeying the d'Alambert equation), 

^(^)(xi) ^e-*"'^(^i'*)'|(^)(xi) , 
^(^)+(x2) ^ e+^'^^(^2,t)|,(e)(x2) , 

A(X2 + ^(X1 -X2),t) ^ A(X2 + ^(X1 -X2),t) - 
dx (X2 + ^(X1 -X2),t) 

9X2 ' 

or, equivalently, this arbitrary phase x(x, t) cancels in (p2[). Such invariance is crucial for 



macroscopic physics since it provides the physical sense of the Wigner function ( 21 ) and 
the conservation of the macroscopic electric current. 

Deriving electron RKE (|l^ we used such gauge transformation while for the weak 
interaction terms we did not. 

Thus, a recipe of the gauge invariance restoration in RKE's dll), ([l|) should be similar 
to the change of arguments in (pl|): we should substitute the generalized (congugate) 
neutrino momentum, 

Q = q + G^^/2cyj(^) (x, t) + A(^) (x, t) , (24) 

for the neutrino gauge invariant Wigner distribution 

/('^)(q,x,t) = /»(Q,x,t) , (25) 
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and the matrix of the generalized momentum 

P = [p + Gi.\/2cy5j('^)(x,t)]5^y - GFV2cA5n^''\x,t){a)^^' , (26) 
for the electron gauge invariant Wigner distribution 

/(^)(p,x,t) = /»(P,x,f) , (27) 

where P = 2~^Tr P. 

The RKE's for gauge invariant (p5|),(p7|) are obtained from ([T5|), (p^) as 

0<5/»(Q,x,t) . <59/»(Q,x,t) ^ dff\Q) _ 
dt ax +^ 9Q 

and 



a^/»(p,x,t) ^ .^ aj/»(p,x,o ^ 

5t 9x 



dP 



where generalized momenta Q, P are given by Eqs. ( p^ and (26) correspondingly; the 

background distributions /q"^ can be changed to /o'^^(q) for neutrinos and fl)^\£p) (19) 
for electrons since we retain the lowest order ~ Gj? in the corresponding weak interaction 
terms. 

Substituting in last RKE's the momenta (|2j) and (p6|), for which we take into ac- 
count both the derivatives (p!6|), ( [l8| ) and the total time derivatives d5j^''^(x, t)/dt = 
a5j(^)(x,t)/at + (nV)(5j('^)(x,t), d5A(^)(x,t)/dt = a(5A(«)(x,t)/5t + (nV)5A('=)(x,t) in 
the neutrino RKE, d5j[f''(x,t)/dt = d5]i^\:x.,i) / dt + (vV)5j^'^^ (x, t) in the electron one, 
and then using the identities 

(nV)<5j(^)(x,t) - V(n5j(^)(x,t)) ^ -[n x V x 5j('=)(x,t)] , 

(nV)(5A(")(x,f) - V(n5A(")(x,t)) = -[n x V x (^A(^)(x,t)] , 

(vV)5j('^)(x,t) - V(v5j('^)(x,t)) ^ -[v X V X 5j(^)(x,t)] , 

(vV)5n('')(x,t)l3(°) - V(vl3(°)jn('')(x,t)) = -(v x V x bW.Jn^'^) (x, f)) , 

we arrive to the final forms of the lepton RKE's. 

Namely, accounting for the definition of the neutrino gauge invariant distribution (25) 
the neutrino RKE (jl^) which is equivalent to the master ^ takes finally in the Vlasov 
approximation the form 

d5&\c^.^,t) 9(^/H(q,x,t) (v^) g/r(q) , 

9t + 5x ^^^^^^ ^""'^'e, dq, + 



^ ' ^e, dq, 



+6F-'{^,t)^:^^^^ = 0, (28) 



where the antisymmetric tensors ^Fj^^ ' (x, t) entering effective Lorentz force terms are 
given by the vector and axial vector currents correspondigly, 

SFjp{^,t)/GFV2cv = ejkiiV x <5j(^)(x, t)); , 

r- (A) (.) d6A^^\x,t) 
V25Fj^\^,t)/GFCA = -V,6A^^\^,t) , 

V26fI^'^ (x, t)/GFCA = e.fcKV x 6A^'^ (x, t))i ; (29) 
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and d/dQ = d/dq from the relation (^). 

First three terms in RKE ( |28| ) were derived in |lO|. In the paper |12| analogous 
vector coupling terms (~ cy for electron neutrinos) were obtained in cold hydrodynamics. 

Analogously for the electron gauge invariant distribution we obtain finally from 
(|20|) (equivalent to the master equation (|l^)) the electron RKE 



9<5/(-)(p,x,i) , 5<5/(^)(p,x,f) 



+::^V(a(p)<^j^^^(x,i)) 



, 



(30) 



where the weak vector term (~ cy) has the Lorentz structure with the tensor components 

in 

_ _5<^jf)(x,t) 



dt 



f-,{weak) 
jk ' 



(31) 



d/dP = d/dp from the relation (p6|); and the axial vector term does not contribute to 
the continuity equation since the last term with the 3-vector component of the four-vector 
o^(p) introduced after (20) , 

me[Ep + me) 

vanishes under integration / (fp{ak{p) /£p)dS^^'^\Ep) /dpj due to the isotropy of the back- 
ground and the odd power of momenta p under the integral. 

Obviously, RKE's (p8|), ( ^0|) obey the lepton current conservation law 



gji"^(x,t) 

dXn 



0, 



(32) 



or the lepton currents are conserved in a magnetized plasma. 



2.4 Dispersion equation for lepton density perturbations in 
isotropic plasma 

Here checking master equations (|9|), (^) against known results for an isotropic plasma in 
the absence of an external magnetic field, Bq = 0, hence neglecting spin waves, we show 
that accounting for the lepton current conservation (^) such observables like the spectra 
of lepton density perturbations do not depend whether we apply initial RKE's for gauge 
non- invariant distribution functions (|^, (|l^), or same equations written in the completed 
forms, 

The situation is similar to the case of standard plasma where the initial RKE for the 
gauge non-invariant Wigner distribution /('^)(p,x, t) is often more suitable for the obtain- 
ing of concrete results than Boltzman equation with the Lorentz force |jl^ . Nevertheless, 
the electric current should be written in the gauge invariant form obeying (^). 
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As well as authors ||^, |To| we consider the electron plasma waves (EPW) driven by the 
intense neutrino flux, neglecting also the transverse wave contribution, i.e. retaining only 
the electric field in the Lorentz force for electrons. 

Using the Fouirer representation in the RKE's (^), (10) one can easily obtain the 
algebraic system for the current perturbations 
5jt)(u;,k) = /d3g(^^/g)5/H(q,k,u;)/(27r)3, 
5ji^\u;,k) = Jd^p{p^/ep)6f'Hp,k,oj)/{27T)\ 



Sn^"^ (w, k) + GfV2cv f<5n(^) (w, k) 



(27r)2 



uj — kn 



d^q nkdf^''\q)/dq^ 



, 



(27r)3 iu-kn 
Sn^"'^ (w, k) [l + Xe(^, k)l + GfV2cv ISn^"^ {to, k) 



d^p kjdf^'''{ep)/dp 
(27r)3 u; - kv 



-^jf (a;,k)fc,- 



d^P Vkdfo'\ep)/dpj 



(27r)2 



w — kv 







(33) 



where n = q/g is the velocity of the massless neutrino and the term Xe{^, k) is connected 
with the longitudinal permittivity, 

, , , , , , 47re2 f d^p k.;d&\e„)/dpi , 

X.(.,k)=e,(-,k)-l = ^/p^^^^4^. (34) 

Note that starting instead of (|9|), ( p^ ) from the corresponding RKE's (|2|), ( pOD derived in 
the previous subsection for the gauge invariant Wigner distribution functions f^'^\Q, x, t), 
f^'^\P,x,t) we obtain the same algebraic system ( |33| ) omitting axial terms (~ ca) for 
Bo = 0. 

Really additional terms coming from the total time derivatives. 



9'5j(^)(x,^) 
dt 

for Q via the second term in (p^), and 

dt 



+ (nV)5j(^)(x,t) 



+ (vV)5j(^)(x,t 



for P via the second term in (|26D , are proportional in the Fourier representation to the 
Cerenkov factors uj — kn, iv — kv correspondingly that exactly cancels with the poles 
within integrands in ( p3D leading to zero contributions in brackets for both equations 
(H): ~ Sj^'^u;,],) J{d^q/i27Tf)dfi''\q)/dqj = 0, and 
~ 6j(-Huj,k) J{d^p/{27rf)dft\ep)/dp, = . 

Now using the symmetrical tensor entering the second line in (|3^ ) (and the analogous 
one in the fourth line) 



d^q nkdf^''>{q)/dq. 
(27r)3 a; - kn 



A5kj + B- 



■ kkkj 
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where the scalars A, B are given by 



A 



B 



(27r)3(a;-kn) r ' ^q^ P ' dqi 



dh. 



(27r)3(w-kn) \ A;2 ' dq, 



dqi 



we obtain the factor multiplying the electron current 6j'j^^ in the second line (and the 



analogous one in the fourth line in (pj^)) 
d^q nkdf^''\q)/dqj 



kj 



(27r)3 
Substituting eqs. 



kk{A + B) 



kkOJ 



dS 



w-kn ' k^ J (27r)3(w-kn) ' 

p^ ), p5| ) we rewrite the system (p^) as 

^n^'^) (c^, k) + GfV2cv \6n^^^ (u, k) - ^kkSji'^ (uj, k) 



9gi 



(35) 



d^q 



ki 



(27r)3(a; - kn) dq. 



Gi.V2cy[(5nM(c^,k) - ^kk6jj^\u;,k) 



k'^Xe{uJ,'ii) 
47re2 



+ 



+ 



<5n(^)(c^,k)(l + xe(c^,k)) = . 



(36) 



Using the lepton current conservation (p2|), A;jt(5j^"'* (w, k) = u;5n^") (w, k), we obtain from 
I) the dispersion equation in the lowest order over the Fermi constant (~ G|^)0, |l^] : 



/S.^k'^Cy 



UJ. 



pe 



^^^""'^V (27r)3(^-kn)^^^^ 



where A^, = 2Gpneon,^o / rrie is the weak parameter introduced in /q*^^ = f^' {q) /riyQ 
is the normalized neutrino background distribution function; = \/47raneo/"i-e is the 
non-relativistic expression for the plasma frequency. 

Let us stress the importance of additional terms in RKE's (l28|), ( ^0|) discussed after 
eq. (|3^) which did not contribute to (|3^) while they become important when we claim 
the lepton current conservation deriving dispersion equation (^7|). We can conclude that 
the master RKE's (|9|), ([To| ) plus the lepton current conservation laws provided by corre- 
sponding eqs. (|2^) , (l30|) lead to the correct issues known in literature, in particular, to 
the dispersion equation (^7|). 

The prediction of the streaming instability driven by the neutrino beam in a dense 
plasma (i.e., outside the neutrinosphere of a supernova) was criticized in [1lO| , 
In this section, however, we do not touch these issues following from eq. (p7[). 
Let us turn to the case of magnetized plasma. 



(37) 



3 Bloch equation in the presence of neutrino beam 

Substituting the linear decomposition (^ into eq. (p!^), then integrating the latter over 
d^p and multiplying by ^jlb we have generalized Bloch equation for the magnetization 
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density perturbations m(x, =| hb \ J d^p6S{p,x,t)/{27r)^ in a polarized NR electron 
gas that is the base of the theory of the electron paramagnetic resonance in the absence 
of weak interactions [pof while in SM with neutrinos such equation takes the form 



9mj(x, t) 
di 



+ 



(m(x,t) -xob(x,t)) X Be 



j dxk 



nip 



(27r)3 dxj 







(38) 



Here xo = ~'^I^'b I d^Pi^f^'^^ / dsp) / {2tt)^ is the static susceptibility of the polarized electron 
gas 0; b(x, t) is the magnetic field perturbation in the total field B(x, t) = Bq + b(x, t); 
"-eo = / d'^ p f^^^Ep) / {2Tr)^ is the density of the electron gas. The pseudotensor T,kj{x, t) = 
fj^B J d^pvk6Sj{p, X, t) /(27r)^ describes spatial inhomogeneity of the magnetization and the 
last new vector term (~ ca) corresponds to the parity violation for the evolution of the 
macroscopic axial vector mj in SM. 



Note that last two lines in the complicated spin RKE (12) are the relativistic cor- 
rections and do not contribute to ( |38| ) in NR plasma. Moreover, we omitted the small 
term — Gi?\/2cyiXo[j^'^H^) ^ ^o]j originated from the third line (without derivatives) in 
(12) when comparing it with the third term in the Bloch equation ( p8| ) stipulated by the 
standard spin precession in QED plasma (arising due to the last term in the second line 
in (PD). 



We neglected also the fifth term in (12) that is proportional to the vector coupling 
~ cy since the background polarization is small in our WKB approximation, S^'^'^^ ~| 
{lJ,BBQ/'y)df^^^^ /dsp f^^^\£p), and this term is much less than the previous one in (12) 
which is proportional to the axial coupling ~ ca retained in (35). This is in agreement 
with estimates of the small polarization < A >= neo/rie ~ 0.01 made in Q (Nunokawa 
et al) even for strong magnetic field in SN where noe is the electron density on the main 
Landau level, He is the total electron density. 

In order to break the chain of hydrodynamical equations for magnetization moments, 
rrij, Sjj, etc, one can consider long- wave spin waves with the wave lengths As that are much 
bigger than the electron Larmour radius, A^ ^ leH- In such case the spin perturbation 
function is given by the first two moments: 



liBSSj{p,x,t) 



df^ 

de.n 



d^p df^^^ 
(27r)3 dep 



1 -1 



X m 



(x,t) + ^Sij(x,i)^ 



that allows us to complete the hydrodynamical system by the equation for the pseudoten- 
sor S,- 



Ft 



+ U,e{ejtinfT.it{yL,t) - euihf T^tj {-x:, tfj + 
-(mj{x,t) - xo&j(x,t)) - 



<v'^>d 



3 

GpCAP-BneO 



nip 



d^q . , a5/(^)(q,x,t) 



(2vr)^ 



r%(g) 



dx 



. 



(39) 



3 



^ The static susceptibility is small in a degenerate electron gas , xo = avp^l^'^'^ <C 1 , in contrast with the 
varying one x(i) (see below eq (pl|)). This is the reason why the static magnetic induction i?o = (1 + 47rxo)-ffo 
and the magnetic field strength iJo practically coincide there. 
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Here Qf, = eBo/rrie is the electron cyclotron frequency; 

< >= [/ d^pdf^'^^/dep]~^ J d^pv'^df^^y dsp is the average of the velocity squared; = 
B0/-B0 is the unit pseudovector and last true tensor term describes the parity violation 
through weak interactions. Note that eq. ( |39[ ) would be important accounting for the 
large spatial dispersion k < v >^ uj outside the region (^) and claims an inclusion of 
exchange interactions (from the ee-scattering exchange Feynman diagram) we omitted 
here. 



4 Neutrino driven streaming instability of spin 

waves 

In this section we derive from generalized Bloch equation ( ^8| ) the dispersion equation in 
magnetized plasma in the presence of neutrino beam. 

Apparently master RKE's for number density distribution functions (P), (|To| ) should 
be consistent with the spin RKE (|l2| ) and its consequence (|38|). Therefore we are check- 
ing below whether more general (P8|), ( |30[) are necessary. As we find below there is no 
difference between the use of these equations for the spin wave propagation in plasma. 

Let us consider for simplicity the case of long-wave perturbations with the spectrum 
a;(k) obeying the inequalities 

— >w»A;<t'> . (40) 

As the mean electron velocity < v > \s small in NR plasma, < f ><C 1, and keeping in 
mind the Maxwell equation written in the Fourier representation b = [k x (5E] /lo one finds 
from the condition k/oj ^< v > that the electric field contribution in the spin RKE (12) 



(~ 5E) is negligible comparing with the magnetic one even for the maximum frequency 
in whole space-like region < k relevant to the Cerenkov resonance with neutrinos, 
UJ = kn. Without electric fields the transversal components of the permeability and 
the susceptibility tensors (i,j = x,y) appearing in the second term of ( |3^ ) are diagonal 
in plasma, ^jj(a;,k) = ^{uj,'k.)5ij, Xiji^-,^) = x{'^-,'^)^ij-, that differs this medium from 
ferromagnets. 

Moreover, the latter inequality in (|4^), uj ^ k < v >, means the high-frequency 
approximation, for which the spatial uniform susceptibility x(^) the permeability 
^{t) = 1 + 47rx(i) can be considered instead of more general ones, x(x, i) and fi(x,t), 
that allows us to neglect the complicated pseudotensor term Sjj in ( |38|) as well as the 
spatial dispersion in the perturbations m±(x,t) = J dt'x{t — t')h±{x,t') and 6±(x, i) = 
/ dt'fi{t — t')h±{x, t') correspondingly where /i±(x, t) = hx{^, t) ± ihyix, t) is the magnetic 
field strength perturbation. 

Under such conditions we find from the linearized spin RKE (|l2|) the susceptibility in 
QED plasma neglecting neutrinos, 

UJ ± ile{l - 47rxo) 

that may be finite at the paramagnetic resonance uj = ^F^^e given by the equation 1 + 
iirxi^^) — 0. Latter follows from the shortened Maxwell equation k x (b — 47rm) = 
when one neglects electric field terms. 
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Then substituting into ( |38D the solution of the neutrino RKE (^) obtained in the same 
Unear approximation (|l]), 

*/<')(,.k,.)=0.^c.^^!^/|^(I-v.„W>.)(p..k..) + 

^ Gfca kkUiiq) df^''\q) f d?Pe , . . 

+ V2 (^-kn) da, J (2^^^'(P-^''^) ' (42) 

we can rewrite the generalized Bloch equation (^8|) in the Fourier representation as 
— iujmj{]i.,uj)+ 2/i5[(m(k, w) — xob(k, w)) x Bo]j H — A\'^ {uj,'k)mi{'k,uj)ikj + 



Svrme \ uj 
x(^eikikk[bi{k,oj) - 47rmz(k, w)] + wJE^) =0 . (43) 

Here the factor is given after eq. (|37|); the vector aI'^^ = A^'^\uj, ICj and the scalar 
= B^'^\uj,'k) depend on the neutrino background distribution /o'^^(q), 

J.)^ f d^q ft\^) / [(kn)^-fc^]nKg) , (kn)n, - fcA 
' 7 (2^)3 g 1^ (u;-kn)2 ^ u;-kn ^ ' ^ ' 

^ / ^ ^(q) ( ii^nr-k'] ] .45. 

Obtaining last term in eq. ( |4^ ) (~ c^cv) we used the exact Maxwell equation (^) when 
we substituted the convection current 6ji(k,u>) = J d^pvi5f^^\p,]i,LL>) in the first line of 
eq. (H). 

Note that the electromagnetic current jjl^^ (x, t) conserves automatically since we can 
neglect weak interactions in the electron RKE ( [To| ) while retaining them in the neutrino 
RKE (P). In turn the neutrino influence the electron spin evolution coming from the 
solution (B^) is not changed if we would substitute the solution of the complete equation 



(28) instead of the initial one (y). 

Really in the last case there appear the two additional terms in (^), 

^ /7T dfo\}x-(e)f IN Gf dfl[\q) 

- GfV2cv — 7; SjI '{uj,k) y=CA — 7; dmk{uj,k) (46) 

oqk v2 OQk 



which do not contribute to the generalized Bloch equation (38). 

One can easily see that in NR plasma with a low electron density noe *C ~ 
2 X 10^^ cm^3 ^]2e general condition of the macroscopic description 27r/k ^ {no)~^^^ 
means that long wave lengths exceeding the Compton one, 27: /k ^ fn'^^ are possible 
only, and due to (|40| ) frequencies obey the inequality lo <k <^ mg. 

Accounting for low oj and k and comparing the term ~ cyCA with the previous one 
in eq. (|4^ ) (~ c^^) one finds that the convection current gives a negligible contribution. 
Note that for f^^,-neutrinos with | cy |<C 1 such correction becomes even less than for fg. 
Omitting also the small term ~ xo we arrive to the shortened form of (^), 

A i \ 

- iujmj{k,uj) + 2/iB[m(k,a;) x Bo]j + ^ " A[' {uj,'k)mi{k,uj)ikj = . (47) 
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Thus, from the generahzed Bloch equation (^) we have derived finally the dispersion 
equation in SM, 

(-^ - ^l) - ^4^)^.1 - ^^^^{A^^'K + A^:'k.) = , (48) 



with the vector A^'^^ as given in eq. (44) 



In the particular case of the neutrino beam, /o'^''(q) = (27r)^5(^)(q — qo), no = qo/'ZOi 
for the transversal neutrino propagation no -L Bq with the beam direction along x-axis. 



no = (1, 0, 0), and the magnetic field Bo = (0, 0, Bo) we obtain from eq. (|48D the resonant 
excitation of spin waves uj = kno + i(5 = f^g + i(5 with the increment 5, 

^2/3 




A;sin0,j2/3 > {V2\ca\ sin^^J 

4 \ go / 



(49) 

where we substituted the scalar product kno = k cos Og^ = k sin 6 cos (p denoting Og^ as 
the angle between the neutrino beam direction and the wave vector k; 9 is the angle of k 
with respect to the magnetic field Bo- Note that we relied in the last inequality on the 
frequency approximation A; > a; ~ assumed above. 

Let us compare this increment with the fastest one in the case of isotropic plasma 
following from the dispersion equation (^) 0], 



V3 f A 

V sm 



6, 



1/3 



2 ' Kqo cos4 



10 . 



One can easily see the advantage of the excitation of collective modes by the intense 
neutrino flux in a polarized electron gas. For a strong magnetic field in a dense plasma 
obeying > ujpe the angular dependence in eq. (^9|) ~ (sin Og^)'^/'^ gives a less suppression 
of the increment for the small angles Og^ < arccos(< v > /c) for which Landau damping 
for growing modes is absent. This is due to the absence of the factor (1 — w^/fe^)^ in 
the new dispersion equation (|48|). It is obvious that such angular dependence would be 
especially dangerous and important for the relativistic plasma case. 

There is the second advantage of spin waves enhanced via the weak axial vector cur- 
rents (~ CA = ^0.5) instead of the case of plasma waves excited via the weak vector 
currents with the small vector coupling in the case of muon and tau neutrinos (choosing 
the lower sign in cy = 2^ it 0.5 where ~ 0.23 is the Weinberg parameter ). This is 
the reason why authors Q considered the case of electron neutrinos only and put for 
them Cy — > 1. Note that during the main neutrino burst in SN all neutrino species are 
produced in the hot SN core via the pair annihilation e^e~ — > z^a^^a; o = e,/i,r. 

On the other hand, there are other arguments against streaming instability driven by 
neutrino beams in SN (e.g. not collimated beam) ||l^ to be relevant for neutrino prop- 
agation in a magnetized medium. Nevertheless, we have just showed that in a polarized 
electron gas the dispersion equations are quite different from the case of the isotropic 
plasma [0, 10, |l^ that stimulates a future exploration of collective plasma phenomena in 



the presence of intense neutrino fluxes. 
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5 Discussion and conclusions 



The above derivation shows that the Bogolyubov method starting from QFT Feynman 
diagrams remains a power and straightforward tool to obtain RKE's. Being supported 
with an additional gauge restoring transformation similar to |17| it leads to the same 
results obtained within Hamiltonian approach of 0] for the case of neutrinos propagating 
in an isotropic plasma. 

In the case of a polarized electron gas the master RKE's derived above by the Bo- 
golyubov method allow to analyse a new phenomenon- spin wave propagation in NR 
plasma enhanced by the neutrino beam. Note that complete RKE's (||), (|0|) are not 
necessary to derive the dispersion equation (^) for spin waves as we showed in (p6|). 

The violation of parity in the SM lepton plasma given by axial vector currents (~ ca) 
leads to the growth of spin eigen modes through their excitation by the intense neutrino 
flux. These spin waves are generally coupled to the magnetosonic ones analogously to the 
case of spin waves in ferromagnets [pH or can transfer their energy to electromagnetic and 
plasma waves at the cross of spectra that finally could lead to the heating of ions and the 
background plasma. 

The possible explanation of the shock revival in SN by different collective mechanisms 
including the neutrino driven streaming instability seems to be very perspective goal for 
future studies in the case of polarized electron gas. 

The application of these mechanisms in the magnetized plasma behind the shock and 
outside the SN neutrinosphere is self-consistent with plasma parameters expected there. 
Really, we do not consider neutrino collisions with matter within this region using the 
collisionless neutrino RKE and may also neglect the electron-ion collision frequency 
fgj comparing with the cyclotron frequency at the paramagnetic resonance lo = Q^. Hence 
collisionless Vlasov approximation should be valid for the spin equation (^) as well as for 
the electron RKE (0) since the Debye number is large, Nd = nQe.r% ^ 1. Indeed, in the 
field Bq = 10^^ Gauss the cyclotron frequency reaches Qe = 1.7x 10^i?o ~ 1.7x 10^^ sec~^. 
This frequency is comparable with the plasma one at the density ngo ~ 10^^ cm""^, 
LOpe = 1.8 X 10^^ sec~^, and turns out to be larger than e.g. the electron-proton collision 
frequency Ugp = 50neo{cm^^) / {Te{K))^/'^ sec^^ ~ 1.6 x 10^'^ sec~^ in the surrounding 
NR plasma with the temperature Te 10^ K. For these parameters the Debye radius 
r£) ~ 10~^ cm corresponds to the large plasma parameter Njj ~ 100 ^ 1. 

Under such conditions for the mean neutrino energy go ~ 10 MeV the increment (p9D 
reaches the maximum value 6 ~ lO^'' sec~^ as in previous optimistic estimates[0, p!o| . 
Note that this means too sharp collimated neutrino beam with the spread of directions no 
for the fixed wave number k not exceeding very small value A(kno)/fie ^ S/Qe ~ 10~^. 
If neutrinos move along radii beyond the neutrinosphere r > Ri, this estimate is too 
optimistic since there is a spread of the angular distribution of neutrino trajectories that 
damages simple model with parallel rays assumed here [10|. 

Nevertheless, we think the simple dispersion equation (48) based on the enhancement 
of pure magnetic field perturbations is the only particular case of the general kinetic 
equations derived here. 

Note that even in this approximation there are some advantages of our model com- 
paring with the isotropic plasma case discussed in g |l|, |l|] as we have just shown in 
previous section (after eq. (|50|)). 

A more general case with the overlap of electromagnetic eigen modes in a polarized 
medium and accounting for the spatial dispersion in such plasma seems to be a more 
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realistic model for a magnetized SN while it is beyond of the scope of the present work . 

We would like to acknowledge discussions with A.S. Volokitin and A.I. Rez. This work 
was supported for V.S. by the RFBR grant 00-02-16271. 
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